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Two crucial issues in the interface studies with the interface specific nonlinear spectroscopy are on 
how to describe the macroscopic dielectric constant of the molecularly thin layer and how to calculate 
the anisotropic two dimensional microscopic local field effect. In the currently accepted models of 
the nonlinear optics, the nonlinear radiation was treated as the result of an infinitesimally thin 
(3JT)' polarization sheet layer, and a three layer model was generally employed. The direct consequence 

of this approach is that an apriori dielectric constant, which still does not have a clear definition, 
has to be assigned to this polarization layer. Because the Second Harmonic Generation (SHG) 
and the Sum-Frequency Generation vibrational Spectroscopy (SFG-VS) have been proven as the 
sensitive probes for interfaces with the submonolayer coverage, the treatment based on the more 
realistic discrete induced dipole model needs to be developed. Here we show that following the 
molecular optics theory approach the SHG, as well as the SFG-VS, radiation from the monolayer 
or submonolayer at an interface can be rigorously treated as the radiation from an induced dipole 
lattice at the interface. In this approach, the introduction of the polarization sheet is no longer 
necessary. Therefore, the ambiguity of the unaccounted dielectric constant of the polarization layer 
is no longer an issue. Moreover, the anisotropic two dimensional microscopic local field factors can be 
, explicitly expressed with the linear polarizability tensors of the interfacial molecules. Based on the 

planewise dipole sum rule in the molecular monolayer, crucial experimental tests of this microscopic 
' treatment with SHG and SFG-VS are discussed. Many puzzles in the literature of surface SHG and 

^ , SFG spectroscopy studies can also be understood or resolved in this framework. This new treatment 

• may provide a solid basis for the quantitative analysis in the surface SHG and SFG studies. 
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I. INTRODUCTION 



The past two and half decades have witnessed tremendous advancement and applications of the in- 
terface specific second order nonlinear optical techniques, mainly the surface second harmonic gen- 
eration (SHG) and sum- frequency generation vibrational spectroscopy (SFG-VS), in molecular interface 
studies^^^^^. 11 ! 12 ! 13 ! 14 ! 15 ! 16 . 17 ^ 8 ! 19 . 20 ! 21 ! 22 ^^ To simply put it, SHG is the second order nonlinear process 
where two photons with the same fundamental frequency (ui) interact with a nonlinear medium simultaneously to 
generate a photon with the second harmonic frequency (2u>). If the two fundamental frequencies are not the same, a 
photon at the sum of these two frequencies can be generated from the so-called SFG process. Because of the sym- 
metry requirement for the second order nonlinear processes, the leading dipolar term of the SHG or SFG processes 
is generally forbidden for the centrosymmetric bulk medium, and thus SHG and SFG become effective probes for the 
interface between the two centrosymmetric bulk phases.— The theoretical foundation and the experimental demon- 
stration of the interfacial selectivity of SHG as well as SFG were pioneered by Shen and his colleagues since the early 
1980'S) 27 i 28 ' 29 i 30 ' 31 i 32 ' 33 ' 34 extending from the original formulation of the light waves at the boundary of nonlinear 
k>( | media by Bloembergen et al. in the early 1960'si 3 ^ 

With SHG and SFG-VS, equilibrium and dynamic behaviors of the molecular interface or film can be di- 
rectly measured from the nonlinear electronic or vibrational spectroscopic response of the interfacial molecular 
moieties i 2 -^ ' 10 ^ 1 ' 24 One particular aspect of the studies in the past few years has been focused on the quantita- 
tive measurement and interpretation of the molecular orientation and vibrational spectra from the SHG and SFG-VS 
measurement on various molecular interfacesi 21 i 25 ' 36 i 37 Recent works also demonstrated that the coherent nature of the 
surface SHG and SFG-VS processes makes them more advantages over the other incoherent spectroscopic techniques 
used in surface studies. The interference of the molecular electric field and the strong polarization dependence in the 
SHG and SFG-VS response can be employed to investigate the detailed interactions and to determine the molecular 
conformation at the molecular interfaces i 38 i 39 ' 40 

However, researchers have long realized that the crucial issues in the quantitative interface studies with the SHG 
and SFG are on how to describe the macroscopic dielectric constant of the molecularly thin layer and how to calculate 
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the anisotropic two dimensional microscopic local field effect ^L^L^^^^L^^L^^. It has been well 
demonstrated that in majority cases the calculation of the molecular orientation is quite sensitive to the values used 
for these two factors.^^^^^ 

In the currently accepted models of the nonlinear optics of interface ; 4 ' i 1 '' 4 ' 55 the nonlinear radiation was treated 
as the result of an infinitcsimally thin polarization sheet layer, and as the starting pint a macroscopic dielectric constant 
was assigned to this thin layer . 4 ' 34 ' 54 ' 55 This model was a natural extension of the original formulation by Bloembergen 
and Pershan, where a nonlinear plane parallel slab with a finite thickness embedded between two linear dielectrics was 
the source of the SHG radiation, and where the Maxwell's equations which satisfy the boundary conditions at the two 
plane interfaces were accurately solved* 3 - 5 " However, according to that treatment, Bloembergen and Pershan concluded 
that the surface dipolar contribution to the SHG signal should be overwhelmed by the quadrupolar radiation in the 
much thicker boundary layer. Later experimental observations demonstrated that the surface dipolar contribution 
could be dominant in many cases. By solving the Maxwell's equations which satisfy the boundary conditions at the 
infinitcsimally thin polarization sheet layer, Heinz and Shen provided the theoretical basis for the interface specific 
SHG, as well as the SFG-VS, in the interface studies. 4 ' 34 ' 54 ' 5 ^ 

Parallel to this macroscopic treatment, Ye and Shen employed a classical induced point-dipole model to include the 
microscopic local field effect on the nonlinear optical properties of adsorbed molecules on a substrate. 56 Shen later 
realized that from the theoretical point of view, the dielectric constant in the macroscopic model is not well defined 
for a monolayer because it is a macroscopic or mesoscopic property. Thus, Shen et al. tried to phcnomcnologically 
interpret this macroscopic dielectric constant as a result of the microscopic local field correction in a monolayer, and 
they also gave the explicit expressions for the macroscopic Fresnel factors and the microscopic local-field factors . 41 i 53 
Zhuang et al. also demonstrated that in a Langmuir monolayer it is satisfactory to treat the whole molecule with 
one unique microscopic local field factor derived from a modified Lorentz model of the interfaced This approach has 
been widely used for interpretation of the SHG and the SFG-VS data sincei 21 ' 25 

However, puzzles and confusions remained in the SHG and SFG-VS practices on how the linear macroscopic 
dielectric constant and the microscopic local field factors, as described by Ye and Shen, should be used in the 
quantitative treatment of the experimental data. In their SHG measurement of the self-assembled-monolayer (SAM) 
at the gold substrate, Eisert et al. carefully compared the results with the macroscopic three layer model and the two 
layer model, as well as the local field corrections. They concluded that using the two-layer model without local-field 
correction gave most satisfactory agreement of the molecular orientation with the results from the NEXAFS and IR 
spectroscopy measurements* 5 ^ This work is obviously not consistent with the above treatment by Shen et al. It was 
further criticized by Roy on the inconsistency that the SAM film is considered anisotropic for SHG but assumed 
to be isotropic in terms of its linear optical properties* 5 - 7 - Roy then outlined a macroscopic phenomenological model 
that treats linear and nonlinear optics of the anisotropic interface layer in a single framework, which nevertheless 
requires a minimum of six unknown parameters of the interfacial optical constants even for an uniaxial monolayer in 
SHG or SFG-VS. 57 ' 58 Thus, Roy concluded that the conventional SHG measurement of phase and intensity may not 
be enough to determine the orientation of SAM. However, even though the anisotropy is undoubtedly significant on 
the nonlinear responses, Roy's criticism might have exaggerated the influence of the anisotropy on the linear optical 
property within the monolayer, which in fact only slightly modifies the reflectivity of the surface as judged from 
ellipsometry measurement ! 26 ' 59 ' 60 ' 61 

There have been a few experimental and empirical theory studies on the SHG data from the interfacial 
monolayer, using different parameters of the macroscopic three layer model plus the microscopic local field 
factors. 43 ' 44 ' 45 ' 46 ' 47 ' 48 ' 49 ' 50 ' 51 ' 62 ' 63 ' 64 In dealing with the macroscopic dielectric constant, generally one of the two 
bulk phase values were used, reducing the macroscopic problem to the same as the so-called two phase modeh 65 ' 66 
However, because the two phase model does not consider the different microscopic properties of the surface molecular 
layer, it can only phenomenologically describe the surface SHG process, and it generally fails to qualitatively explain 
the observed differences in the SHG measurement. On the other hand, in dealing with the microscopic local field 
factors, some works followed the classical dipole model as presented by Ye and Shen, and the rest followed the simple 
Lorentz-Lorenz local field expression using the bulk polarizabilities. 63 ' 64 

Another approach completely neglects the treatment of the microscopic local field effect, and treated the molecular 
film with a macroscopic three layer modeh 42 ' 52 ' 58 ' 67 ' 68 ' 69 ' 70 ' 71 In SFG-VS, some considered the film anisotropic and its 
optical constant is determined using the Clausius-Mossotti relationship with the ellipsometry data of the film* 7 ^ This 
is obviously rooted in the practices of the ellipsometry studies, where it is generally believed that the macroscopic 
optical model can be valid down to the monolayer leveL 59 ' 60 In SHG, some considered the film isotropic and the optical 
constant can be obtained from the Kramers-Kronig dispersion relationship measured from the bulk or thin film UV- 
Visible absorption spectrai 42 ' 69 ' 70 When the molecule is considered with uniaxial symmetry, the ratio between the 
dielectric constants of the fundamental and the SH frequency can be directly determined from the intensity ratio in the 
polarization measurements. 37 ' 42 There are also efforts trying to show that the dielectric constant of the extremely thin 
monolayer is just the simple arithmetic average of the two bulk dielectric constants. 52 ' 58 ' 67 ' 68 Event though the results 
from these approaches seemed reasonable, neglecting the anisotropic microscopic local field effect in the molecular 
monolayer is not likely to be physically sound. 

All the above approaches gave the similar trend in the variation of molecular orientation. However, the value of the 
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oricntational angle depends quite significantly on the different parameters used. Besides all the different approaches 
in the SHG and SFG-VS data analysis, the experimental studies have shown that the initial treatment by Heinz and 
Shen with the infinitesimally thin polarization sheet layer is generally correct. Researchers in the field also agree that 
the definition of the dielectric constant of the monolayer or even submonolayer film is unclear and meaningless in the 
macroscopic optics, because when the film is as thin as one molecular layer the macroscopic reflection and transmission 
coefficients simply converge to the bare interface limit . 58 ' 59 ! 60 Therefore, even though the treatment provided by Shen 
and coworkers is generally valid, there are still issues need to be clarified. One problem in the existing model is that 
the macroscopic and microscopic linear optical properties of the nonlinear molecular layer are not totally compatible. 
In practice we also need a more detailed microscopic theory to help understand the detailed SHG or SFG-VS responses 
from different molecular interface or film. 

Aside from the above, a microscopic theory of molecular crystal surface SHG was proposed by Munn in the early 
1990's^ In this formulation, since the surface SHG response is treated microscopically as a sum of responses from 
successive surface layers, the introduction of a dielectric constant of the surface layer was not required, and this 
incidentally yielded a microscopic expression of the dielectric constant of the surface layer. Subsequent studies by 
Munn and coworkers used the planewise sum rules in the molecular crystal dielectric theor y 73 1 74 1 75 1 76 1 77 1 78 to simulate 
the linear and nonlinear optical responses of the model Langmuir- Blodgett films 49,79,80,81,82 In thege treatmentj 
Munn et al. concluded that improper treatment of the local fields can result in significant errors in determination 
of the molecular angle from the SHG datai 79 ' 80 They also showed that in order not to overestimate the microscopic 
local field factors for the closely packed monolayer films, the monolayer itself had to be segmented into several 
layers ! 49 ' 82 In the microscopic theory, the microscopic local field factors depend on the orientation and distribution 
of the interfacial dipoles, and the microscopic local field factors are needed to determine the dipole orientation from 
the SHG data. Therefore, Pannhuis and Munn suggested that a self-consistent approach needs to be employed to 
solve the problem^ These microscopic theory and simulations certainly provided new insights on the treatment of 
the SHG from the molecular monolayer as well as multilayers. These insights may contain answers to the questions 
raised above. However, the implications of these works are yet to be picked up by the field. 

Because the SHG and the SFG-VS have been proven as the sensitive probes for interfaces with the submonolayer 
coverage^&iii the treatment based on the more realistic discrete induced dipole model needs to be developed. 
According to Born and Wolf, the molecular optics theory can directly connect the macroscopic optical phenomena 
to the molecular properties, and can provide deeper physical insight into electromagnetic interaction problems than 
does the rather formal approach based on Maxwell's phenomenological equations! 83 ' 84 From the above analysis, we 
realized that different from the infinitesimally thin polarization sheet layer model treated on the basis of the Maxwell 
theory, a molecular optics treatment can be developed to describe the coherent SHG or SFG-VS radiation from a 
discrete lattice in between the two isotropic phases. We shall show in this report that even though there is no general 
methods to solve the integro-differential equations of molecular optics, the summation, or integration, of the radiation 
in the far field from a discrete induced dipole lattice can be rigorously solved with the application of the principle of 
the stationary phased 4 - 

The discrete induced dipole lattice model is a more realistic description of the molecular interface than the in- 
finitesimally thin polarization sheet layer model. 85 Because the SHG and the SFG-VS processes are known with 
interface selectivity, the whole problem of the integro-differential equations is greatly simplified since the summation 
over the second harmonic or sum-frequency radiation fields from the discrete interface induced dipoles needs to be 
calculated. According to the Ewald-Oseen extinction theorem in the molecular optics, there exists only the incident 
field in vacuum and the dipolar microscopic radiation field in vacuum emitting from each induced dipole i 35 ' 83 ' 84 Since 
the whole problem is treated microscopically as radiation from the discrete induced dipoles, and the total radiation 
is calculated through summation, or integration, over the whole dipole lattice, there is no need to use the Maxwell 
equations together with the boundary conditions as in the infinitesimally thin polarization sheet layer model. This 
works because as Lalor and Wolf pointed out that the Ewald-Oseen extinction theorem just plays the role of the 
boundary condition^ Therefore, there is no need to use the boundary conditions, where a macroscopic dielectric 
constant has to be defined in order to apply the boundary conditions for a finite volume across the boundary area 
which includes the molecular monolayer. Therefore, the macroscopic dielectric constant for the molecular monolayer, 
which has difficulties to be defined in the previous treatment, is just a nuisance parameter in this microscopic, and 
it is finally able to be gotten rid of. There was a section in their classic paper where Bloembergen and Pershan^ 
showed that the approach using the integral equation based on the Ewald-Oseen extinction theorem reached exactly 
the same results for the nonlinear response from the nonlinear plate parallel slabj 35 ' 86 This certainly helps to justify 
the microscopic molecular optics approach to be employed. 

In the next few sections, we shall show that the SHG radiation in the far field from the discrete dipole lattice is in 
the same form as from the infinitesimally thin polarization sheet layer model. Moreover, in this discrete induced dipole 
model, the introduction of the polarization sheet is no longer necessary, therefore, the ambiguity of the unaccounted 
dielectric constant of the polarization layer is no longer an issue. Incidentally, the anisotropic two dimensional 
microscopic local field factors can be explicitly expressed with the linear polarizability tensors of the interfacial 
molecules. Based on the planewise dipole sum rule in the molecular monolayer, crucial experimental tests of this 
treatment with the SHG and SFG-VS experiments are going to be discussed. We shall also discuss the puzzles in 
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the literature of the surface SHG and SFG spectroscopy studies as discussed above. This microscopic treatment can 
provide a solid basis for the quantitative analysis of the surface SHG and SFG studies. 



II. THE INDUCED DIPOLE OF THE MONOLAYER WITH THE DISCRETE POINT-DIPOLE 

LATTICE MODEL 

Because the SHG and the SFG-VS have been proven as the sensitive probes for interfaces with the submonolayer 
coverage^&i£ the treatment based on the more realistic discrete induced dipole model needs to be developed. Here 
we derive the expressions of the linear and nonlinear induced dipole of the monolayer with the discrete point-dipole 
lattice model. Some aspects of the discussions below can be found in the earlier paper by Ye and Shen^ 

In the classical molecular optics, the response of the medium to the incident field is described by means of the 
electric-dipole moments that are induced in the molecules of the medium under the action of the incident fieldi 83 i 84 
When an optical field at a frequency u> is incident on a medium, it creates induced dipole at the incident frequency and 
its higher harmonics in that medium through the total field that each dipole or molecule experiences. This total field 
is called the local field E\ oc . The higher harmonic induced dipole can be viewed as the results of multiple interactions 
with the local optical field. Therefore, they are only strong enough for detection when the optical field is intense 
enough. The induced dipole at the incident frequency is the source of the radiation in the linear processes, while the 
others are the source responsible for the radiation in the nonlinear processes. Simply to put it, one has 



Nonlinear optics has been an extensively studied field since the invention of the first laser in the early 1960's i 26 i 
Here we only discuss the linear and the second harmonic process in an optical medium. Considering the fact that both 
fields in the linear frequency and the resulted second harmonic can contribute to the induced dipole at the second 
harmonic frequency, one has 



Here ctij is the linear polarizability tensor and flijk is the second order nonlinear polarizability tensor of the each 
molecule in the laboratory coordinates X(x, y, z) with the index ijk as either one of the three laboratory Cartesian 
coordinates. The first term in the Eq[2]is the linear induced dipole, the second term is the induced dipole at the second 
harmonic frequency (2w)induced by the field with the fundamental frequency u through the second harmonic process, 
and the third term is the induced dipole at 2lu induced by the field with 2uj through a linear process. The field with 2lo 
in the third term is the result of the second term. The molecular polarizability tensors in the laboratory coordination 
X(x,y, z) can be projected from the molecular polarizability tensors in the molecular coordination A'(a,6, c). 

The above is general for any dielectric medium. In the problem we describe below, only the induced dipoles at the 
interface with the fundamental and second harmonic frequency are considered. The case for the sum-frequency works 
similarly following the same approach. For the isotropic molecular monolayer, there are two independent linear optical 
polarizability tensors in the laboratory coordination X(x, y, z), a xx — ot yy , and a zz ; and there are three non- vanishing 
independent f3 elements, /3 XXZ — /3 yyz — (3 XZX — (3 yzy , (3 ZXX — (3 zyy and f3 zzz i 25 ^ 53 Now the induce dipole moment of 
each molecule can be written explicitly as 



r* rflinear , -fnon — linear 

^induced — [^induced ^induced 



(1) 



Mi = Mi +Mi 




(2) 



5 





^xx loc,x 






^y 


®"yy-^loc,y 








^zz loc ,z 








= o?*E} w - 

W IX loc,x 


\~ fixxzElocxEfoc,* ' 


t" PxzxE loc z Ei oc x 




= a 2uJ E^ - 
yy ivc^y 




~ ftyzyEi oc z Ei oc y 




- rv 2uJ 

zz loc.z 


V @zxxEi oc x E loc x - 


t~ PzyyEi oc y E loc y 



PzzzEf ocz Ef ocz (3) 

Now, in order to calculate these induced dipoles of the monolayer, the knowledge of the local fields Ef oc and Ef^ c 
at the interface needs to be described. 

When an optical field is incident on a molecular monolayer at the interface between two bulk phases with dielectric 
constant e\ = n\ and £2 = n|j the induced dipole in the monolayer and the image induced dipole in the substrate 
all add to the local field at each individual molecule within the monolayer. Therefore the total local field is different 
from the applied electric field. Then ; 56 ' 88 ' 89 



Eioc (r) = E{f) + E dip (r) + E dlpJ (f) (4) 

Here, let r be the position vector within the monolayer from the center of the molecule at the origin position. Then 
Eioc(r) is the local field at r, and it is the sum of three contributions. Here E(r) is the applied field, Edi P {r) is the 
field acted at this point from all the induced dipoles except the one at r, and Edi P j(r) is the field created by all the 
image induced dipoles at r. 

Here we first discuss the first term E(f). In this case, a field of plane wave at frequency to is incident at an angle 
of Oii upon the interface from one bulk phase and reflected from the interface. Because the applied field is not with 
frequency 2uj, one always have E{f,2uj) = 0. For E(r,oj), because of the existence of the interface, the applied 
electric field is the superposition of the incident field plus the field reflected from the interface ! 56 ' 88 ' 89 According to 
the Fresnel formulae^ the total applied field E(f,u>) (E x ,E y ,E z ) can be determined from the incident electric field 
E (Eo, x ,E 0t y,E 0tZ )=E e as followings: 
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ldcosn ' j E z>0 = L' zz E zfl (5) 

This may be simplified as E = [L' • e]Eo. Here, fli and fit are the incident angle and the refraction angle at the 
interface. k\ is the wave vector in first bulk phase, d is the distance of the dipole from the interface. Generally for the 
monolayer d <§; A, e l2k i dcosn i ~ 1. In this case, the expression of the Fresnel factor L' u is simplified into La, which 
will be expressed in the Section HTT1 

The expression of E(r) in the Eq|5]is independent from the structure of the interfacial monolayer for any interface 
between two isotropic bulk phases. However, the expressions of the other two terms in the EqU depend on the model 
of the monolayer structure. In order to evaluate them, here a square lattice model is used to represent the molecular 
monolayer. Other kinds of lattice models, such as the hexagonal lattice model, can also be employed and they shall 
generate similar results . 90 ! 91 Here we only discuss the case using the square lattice model. 



Bagchi et a l 88 ' 89 and Ye et al^ used the classical discrete square point-dipole model to discuss the linear and 
nonlinear optical responses of the molecular monolayer, respectively. Such a square lattice model is illustrated in 
Figffl The dipoles at the interface forms a two-dimensional square lattice with a lattice constant a. The two bulk 
phases occupy the two semi-infinite space of z > and z < respectively, and they are characterized by the bulk 
dielectric constant t\ and £2 (the substrate), respectively. The tilt angle of the dipole with the surface normal is 8 
and its azimuthal orientation is randomly distributed in the plane of the interface. The distance of the point dipole 
to the lower substrate surface is d. The coordinates of each point dipole is then R m ,n = (ma,na,d), where m and n 
are integer indexes along the x and y directions from the origin O(0, 0, 0) in the plane of the monolayer, respectively. 

The contribution to the dipolar field Edipif) at an adsorbed dipole comes from the neighboring induced dipoles 
within the monolayer^ 9 - As in previous treatments ; 56 ! 88 ' 89 the near field approximation of the local field is employed. 



Substrate 



FIG. 1: Illustration of an infinite lattice of point induced dipoles on a substrate interface, a is the distance between the nearest 
molecules, or the lattice constant; d is the distance from the center of the point-dipole to the surface of the substrate; 9 is the 
tilt angle of the point-dipoles. 



In the near field approximation, the magnetic part of the electromagnetic field is ignored, and the retardation effects 
are discardedj 77 i 88 ' 89 This assumption is justified because the wavelength A of the incident electric field is much larger 
than the lattice constant a. Following the works by Bagchi, the electric field at the lattice point i?oo(0,0,(i) which is 
created by the rest of the induced-dipoles in the lattice plane i o 88 ' 89 
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The CGS unit system is used in this paper, as in many textbooks. Here //$ is the components of the induced dipole 
fl. x, y, and z are the unit vectors in the laboratory coordinates. The symbol prime ' denotes that the summation 
does not include the point of origin Rqq. The constant £o was evaluated by Topping in 1927^ 
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In the EqO only the induced dipole in the same plane of the interface lattice is considered. To be more rigorous, 
molecular layers other than the interface lattice plane also need to be considered. These layers may be the molecules 
or atoms of the two bulk phases, whose polarizabilities differ from those of the dipoles in the interfacial layers. 
Phcnomenologically, let's define a vertical lattice index I and lattice constant d, with I = for the plane of the 
monolayer rather than the plane with the origin O(0, 0, 0). Then the total E dip (r) = + E*™ M ° , with E*° p °< 1 ^ 
as the induced dipole interaction contribution to the Rq Q = Rqo origin point from the induced dipoles other than the 
interface lattice plane. 
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Here a/ is the lattice constant of the I layer. The expression here is similar to the expressions for the image dipole 



term by Bagchi and Ye et al 
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Here £q can be directly derived from the EqJS] 
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To evaluate how quickly this summation converges with increase of I, it can be converted into a more rapidly 
convergent seriesi 88 i 89 
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The EqQIj] indicates that £q falls off exponentially as the I increases. Generally, when d ~ a/2, even the contribution 
of the immediate neighboring layer, i.e. I = ±1 can be neglected. This is the basis for the planewise dipole sum rule in 
the molecular crystal theory, which shall be discussed in the Section HV I 73 ' 74 ' 75 ' 76 ' 77 ' 78 The planewise sum rule in the 
molecular crystal theory concluded that the contribution of the induced dipoles from neighboring layers contribute 
insignificantly to the local field of the induced dipoles in the layer under consideration, comparing to that from the 
induced dipoles in the same layer. 

The image dipole contribution is just a special case in the EqJ5] Bagchi and Ye et al. discussed and evaluated the 
image dipole contributions previousl y 56 ' 88 ! 89 The image dipoles is located in the substrate at R' mn — (ma,nb,—d), 
i.e. I = —2 in the general case above, and the image dipole is defined a o 56 ' 88 ' 89 
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Thus, the image dipole contribution can also be calculated as the same procedure as in the Eq|8 ] 56 ' 88 i 89 Ye and Shen 
concluded that for the typical case on a metal surface, where (e2 — ex)/ (e2 + ei) ~ 1, when d > 2.5A with a = 5A, the 
image induced dipole contribution is negligiblei 5 - 6 - Therefore, it is generally accepted that for the dielectric substrate, 
where (e2 — ei)/( £ 2 + £i) <C 1, this image induced dipole term needs not to be considered ! 43 ' 44 ' 45 ' 46 ' 47 ' 48 ' 49 ' 50 ' 51 ' 81 
Therefore, we shall neglect the image induced dipole term in the followed treatment. 

It has to be noted that both the Edi P {f) and the Edi P j(r) terms have the same expressions for the fundamental 
and second harmonic frequencies. 

With all the three terms known from above, put the EqHJ the EqJS] the EqJS] into the EqJ31 one has 
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Here the Ef is from the Eq|5j and la is the microscopic local field factors for the isotropic monolayer using the 
discrete induced dipole model, la is explicitly derived and is defined as belowi^ 6 - 
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(13) 



Here u>i = lo or 2u>, and a^/ is the linear polarizability tensor which can be calculated from the molecular polariz- 
ability values in the molecular coordinates frame A' (a, 6, c). The summation term in the Eq|13lis the dipole interaction 
from layers other than the interface layer, which are generally negligible. However, in the Section ITVl the form shall be 
used in the discussion of the segmentation of the chain molecules in the molecular monolayer. The £o = —9.0336 is for 
the square lattice model as derived by Topping^ Since £o is negative, l xx = l yy is generally larger than unity, while l zz 
is smaller than unity. However, when a is small, |a^j.^o/2a 3 | > 1 can happen, this would result in negative l xx = l yy 
values even under the normal dispersion condition for the bulk material. This is unique for the two-dimensional cases. 

The molecular polarizability at the experimental coordinates can be expressed as from the coordinates transi- 
tion relation: = Ru> Rjj>cti>j> , with R\^\> as the element of the rotational transformation matrix from 

i' ,j'= a,b,c 

the molecular coordination A' (a, b, c) to the laboratory coordination A(x, y, z). 21 ' 92 For the ensemble of rotationally 
isotropic monolayer in the x-y plane, each individual molecule possesses the same homogeneous distribution in the 
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azimuthal orientation and the twist orientation. Then, one has 



i • (l + cos 2 6) ■ a aa + | • (l + cos 2 0) -a bb + Y sin 2 I 
| • (1 + cos 2 6) ■ a aa + ^ • (l + cos 2 0) • ahb + ^ • sin 2 ( 



= i • sin 2 • a aa + ^ • sin 2 9 ■ a bb + cos 2 9 ■ a cc (14) 

Here we did not consider the case when the tilt angle 9 is with a distribution. In that case, the EqJU need to be 
put into the form with ensemble average over the distribution of 9. In doing calculations, the polarizability of the 
molecule or molecular groups can be obtained from many compiled sources) 93 ! 94 ' 95 or from direct quantum mechanics 
calculations. 

The local field corrected linear and nonlinear polarizabilities in the Eq |12l are the source of the linear or nonlinear 
radiations from the monolayer at the interface. Their detail expressions will not appear in the next section when we 
try to calculate the total radiation from the monolayer. The detailed expression in the Eq[l2l the Eq[l3l and the 
Eq[14|can be used for direct calculation of the local field factors and the local field corrected polarizabilities in general. 

The expression in the EqfT51 can be slightly different if a non-square lattice model is assumed. For other geometries, 
for example, the hexagonal or cqui-triangular geometry, the calculation of different geometries can be put forward 
according to Topping's treatment^ 1 - 

The Eq|13l clearly indicates that when the distance a between the neighboring molecules in the monolayer becomes 
large, i.e. the submonolayer case, the la value is approaching unity very quickly. This actually defines the meaning 
of the word 'local', by considering how rapid the value £o is converged when doing summation of the dipoles in the 
lattice plane over m and n. According to Topping and Philpott's calculations, the convergence is generally reached in 
the 4th digit before m,n < 10i 77 i 91 For a lattice constant of a = 5A, this means the local field calculation converges 
within 5nm. Since when a become large, the local filed factors in the Eq|13l decays to unity rapidly, therefore, the 
local field can be viewed as localized interactions with few nanometers. The planewise sum rule also restricts the local 
filed effect within only a few Angstroms from the dipole under consideration. These facts provide definition of the 
actual range of the local field effects. 

By introducing the microscopic local field factors, the calculation of the response from an ensemble of induced 
dipoles to an electromagnetic field is reduced to a calculation of the response of the isolated molecule interacting with 
the local field. 



III. SECOND HARMONIC RADIATION FROM THE DISCRETE INTERFACE INDUCED DIPOLES 

IN THE FAR FIELD 

Here we present the calculate of the linear and second harmonic radiation in the far field according to the two- 
dimensional lattice of the local field corrected induced dipoles at an interface between the two isotropic bulk phases. 
The calculation using the square lattice model can also be performed with hexagonal or other lattice models. We 
shall see that since the radiation in the far field is additive for the radiating dipoles, ensemble average treatment of 
the radiating dipole can be straightforwardly implemented. The following derivation uses the CGS convention, unless 
specified otherwise. 

In this problem, the radiation from the point induced dipole pi at the linear and second harmonic frequencies is 
in the upper phase (ei) in Fig[TJ Here the point of origin 0(0,0,0) is at the center of the dipole, different from the 
definition in the previous section. This is for the convenience of calculation in this section. Now, the general field 
radiated from an electric dipole /i along the z direction, i.e. fiz, in the E\ phase is^ 







£l 


E 


Ml ( 




£l 


E v 


= 



{(far) 3 (far) 2 



£i 1 (far) 3 (far) 2 (far) 



cosfle < (* l, *^ fart ) 

a i(ki-?-u>t) 



(15) 



Here fi is the modulus of the induced-dipole, fa is the wave vector in phase 1, r is the distance from the point of 
the dipole to the point of detection at M(r, 9, tp) with 9 as the tilt angle and ip as the azimuthal angle in the polar 
coordinates. If we only consider the far field, which is generally in the experimental measurements except for the near 
field studies, all the high order 1/r terms can be neglected. Therefore, only the 1/r term of the Eg in the Eq [TBI need 
to be considered. 



FIG. 2: Polar coordinates with an electric dipole with moment vector along the z-axis. /i is the modulus of the induced-dipole, 
r is the distance from the point of the dipole to the point of detection, 9 is the tilt angle in the polar coordinates, tp is the 
azimuthal angle. 



In the following, the total macroscopic radiation field at a space point R(0, 0, i?o) m the reflection direction can 
be calculated from the summation of the field of the radiation directly from the individual induced dipole at the 
interface and the radiated field reflected from the interface. In doing so, the amplitude and phase of the radiated field 
at R(0, 0, Rq) need to be expressed into the function of the components of individual induced dipole at the interface, 
before the summation can be carried out over them. 



The connection between the induced-dipole components fi x , [i y and ji z and their radiation field components E x , E y 
and E z at R(0, 0, Rq) is illustrated in Figure[3] The total phase at R(0, 0, i?o) can be calculated by choosing the point 
of origin O(0, 0, 0) as the reference point. Therefore the phase difference of the radiation field at R(0, 0, Rq) between 
the dipole at the point r(ma, na, 0) and the fixed point of origin O(0, 0, 0) can be calculated from the difference of the 
two distances RO and Rr, i.e., r(m, n) — R , with 

r(m,n) = \J (ma) 2 + (na) 2 + (Ro) 2 ■ (16) 

Now, the phase difference of the incident field at O and r on the dipole plane also need to be calculated. When a 
plane wave is incident on the monolayer at the incident angle fij, remembering that the incident plane is xz, and m is 
along the x direction, this phase difference is kf 1 (ma sin f2j). For the case of second harmonic generation, this phase 
difference involves two incoming waves, therefore the total phase difference is 2k\ n (ma sinfJi). 

Therefore, the total phase at the space point R(0, 0, Rq) is the sum of all these phase differences. If we define this 
total phase of the radiation from the induced dipole at r(ma, na, 0) at R(0, 0, Ro), as k\f(m, n), one has 

f(m,n) = J (ma) 2 + (na) 2 + (R ) 2 + ^-(masinfl,). (17) 

k\ 

Here when the radiation frequency is the same as the incoming frequency ui, then rj = 1, and k" 1 — k\, while when 
the radiation frequency is 2uj, then n — 2, and k\ n ^ k\. More generally, for the case of the sum frequency generation 
(lus = lu\ + LU2), the expression of the second term in f(m,n) becomes 

-[k^^sin^ 1 +k\ n (L0 2 )sin{lf 2 ] (18) 



fci (lo 3 ) 

With the phase known, the radiation field at R(0, 0, Rq) generated from the dipole pi(pL x , /i y , fi z ) at r(ma, na, 0) can 
be calculated separately for each according to the Eq|15l with proper projection. 

Then the field components E x , E y ,E z at i?(0,0, Rq) generated by jj, x at r(ma,na,0) are 
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ju x (m,n,Q) fj. y (m, n, 0) ju z (m,n,0) 




FIG. 3: Calculation of the field vector (E x ,E y ,E z ), which is created by fj, x (m,n,0), fj, y (m,n,0), fi z (m,n,0) respectively. 6 
being the tilt angle of the dipole moment vector, tp being the azimuthal angle of r which is the vector from the point of (m,n,0) 
to the point (0,0,J?o), as illustrated in Fig[2] 



fcf sin 2 9 

p ' p ' £ir(m,n) 

eir{m,n) 

E^ z = cos g sin y = - fc? Sln ^ C ° S g f D * ^ e ^/( m .n)- M t) ( 19) 



with 



2 



. (m) 2 + Rl -ma 

sin = — = , cos = — 

Vfrna) 2 + (m) 2 + i? 2 . ^(ma) 2 + (na) 2 + i?„ 

Ro -na 
simp — — = , costp= — = (20) 

Similarly the field components -E^, E y ,E z at i?(0,0, i?o) generated by fi y at r(ma,na,0) are 

r-i rn /i fc 2 sin COS COS (yO f , mWwtl 

" eir(m,nj 
eir(m,n) 

E^ z = ^ooBflain V =-M^^^/^-)-^ (21) 



with 



. ^/ (ma) 2 + i? 2 , 

smfl = — — , costf 



2 



\/(m(i) 2 + (na) 2 + i? 2 , ' ^(ma) 2 + (na) 2 + R\ 

Ro -ma 
simp — — == , cosp — — == 22) 

V(ma) 2 + i? 2 ^ VM'+^o 

Then the field components E x , E y ,E z at i?(0,0, Rq) generated by fi z at r(ma,na,0) are 
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n fc? sin 9 cos 9 cos ip i(k H „i_ tl . t i 

= £^,ecos0cos^ = — i ? r — - n z e^ f{m > n > ut > 

eir(m, n) 

E^ y = ^„ fl cos fl Bin y = - * ? ^ ^ f n g ^ /("■")-*) 

*W - -^, e sing= ^!i /ize i(fa/(^,n)-ct) (23) 

emm, n) 



with 



. V( ma ) 2 + ( na ) 2 /) 

sm — — , cos — — 



2 



^(ma) 2 + (no) 2 + R% \/(ma) 2 + (na) 2 + Rq 

—na —ma 
sm tp = — -^^^=^^^= , cos <p = —j^^^=^^^= (24) 
y (ma) 2 + (na) 2 y (ma) 2 + (na) 2 

Now the total radiation field at i?(0,0, i?o) generated by the whole monolayer at the interface is the summation 
over the whole lattice plane m(-oo,+oo) and n (— oo,+co). 

oo oo 

E x = ^ E (-^W^ + E li y ,x + E ^,x) 
77i — — oo n=— OO 
OO oo 

m— — OG 71 — — oo 

oo oo 

E E + + ( 25 ) 



m — — oo n— — oo 



The summation in the Eql25l for the total radiation field is in principle the discrete form of the integro-differential 
equation in the molecular opticsi 83 i 84 The asymptotic result of this summation can be evaluated by employing the 
condition of the far field. Because Rq is very large as compared with the lattice constant a and the wavelength A, the 
discrete summations in the Eq |25l is asymptotic to a continuous integration over the variable set (x,y) = (ma,na). 
In addition, because Rq is large and so the phase factor k\f(m : n) is large, the exponential factors in the Eq [ll)l |2~T1 
and 1231 oscillate very rapidly and change their signs many times as the point r(ma, na, 0) explores the domain of 
integration-^ Under these conditions, an asymptotic value of the total electric field (E x , E y , E z ) at point i?(0, 0, Ro) 
can be obtained with the application of the following formula, which are derived from the principle of stationary 
phase. The detail of the stationary phase method can be found in classic textbook by Born and Wolfi^ 

Here by defining the variables x = ma and y = na with the lattice constant a as a very small quantity compared to 
the infinity size of the lattice plane, each term in the summation as shown in the Eq 1251 becomes the following integral 
and this integral is then asymptotic to the summation of the stationary phase terms as below. 



4 I" r 9(x,yy klfix - v) dxdy=^Yl I = gfo.frOe*^^ (26) 
a 2 J~oo J-oo ha 2 ^ ^a^.-^l 

Here g(x,y) is the pre-exponential factors in the Eq ll9l [2T1 and [23l (xj,yj) is the points in the whole integration 
domain at which / is stationary, i.e. 

The definition of the other terms are the followings. 



a 



-(3L-*-(S0„--(SL 



and 



-1 if cxjflj > 7 2 , ctj > 



-1 if arfj > 4 , ctj < , (29) 
-i if ajfa < y| . 
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Using the stationary phase condition in the Eq |27l and the expression of / in the Eq |17[ a single stationary point 
(xi,y\) is obtained, and 

x\ = — i? tanSl , y\ = . (30) 

Here Q = fi, is the direct results of the stationary phase condition for the linear radiation, and for the second 
harmonic radiation SI has to satisfy the relationship fci(2a;)sinfi = 2k\ n (uj) sin SI;. Similarly, one can also show 
explicitly that for the sum frequency radiation, SI has to satisfy the relationship ki(u)3) sinSl = fc™(wi) sinSl" 1 + 
fe| n (o>2) sin SI" 2 . These are just the condition for the general law of reflection in linear and nonlinear optics, as given 
by Bloembergen and Pershan, as derived from the Maxwell boundary conditions. 3 — From these results we shall show 
later that at the far field space point R(0, 0, Rq), the radiation field is in the direction of SI, which satisfies the general 
law of reflection in linear and nonlinear optics. It is interesting that these conditions came incidentally as a solution 
from the molecular optics treatment. This fact further illustrates the general equivalency of the macroscopic Maxwell 
equations with boundary conditions and the molecular optics treatment] 83 ' 84 

Now with the stationary phase point value, we have 



f(xi,yi) = i? cosSl 
o-i Ro 



v/jaiA -7i 2 | cos2 n 



(31) 



Here <j\ = 1 because a\ > and — j 2 > 0. Now put all these values into the expressions in the Eqf5H we 
have the radiation field components at the space point i?(0,0, Rq) that are directly radiated from the whole induced 
dipole lattice as the followings. 

E* 1 = 0*x cos - ii z sin SI) e t{kl R ° cos 



a 2 £\ cos SI 



= ^fcif-^sinO + / i z ^) e ^«ocos^ t ) (32) 
crei \ cosSz / 

In order to calculate the total radiated field at R(0, 0, Rq), besides the direct radiation from the induced dipoles at 
the interface, the radiation in the forward direction which is reflected from the interface also need to be evaluated. 
The above procedures can be repeated to calculate the radiated field at the image point Ri(0, 0, —Ro) in the isotropic 
phase E\ from the whole induced dipole lattice. We have 

El 1 = — — iki(n x coaQ + ^ z smn)e^ klRocosQ - u ^ 

E T X _ 2?F q Zl ttH p i(feifiocosn-q>t) 



V 



a 2 E\ cos SI 



E^ = ^-ik x ( n x sinSl + ^^-^) e'^^-n-.t) (33) 



a z Ei \ cos it 

Because the choice of the space point R(Q,0,Rq) and Ri(0, 0, — Rq) is rather arbitrary, and the radiation field 
components in the Eq |32l and the Eq |33l are both directional and the amplitudes of both plane waves are independent 
from the value of Rq, the coherence of the radiation from the induced dipole lattice is conserved at the far field. 

According to the Fresnel formulae^ 4 - the amplitudes of the reflected electric field components are 

E R2 = _ r e t, = »icos^t - na coBfli £Tl 

x p x n 2 cos Hi + m cos Sl t x 

R2 _ T _ ni cos fli - n 2 cos Cl t T 
& v — r s& v — tz — ; FT A, 

y y m COS \li + 712 COS S2( " 

^ a cos SI, — ni cos SI* „ T 

E?* = rpET 1 = — -E^ 1 (34) 

ri2 cos ili + n\ cos lit 

The phase difference between the direct radiation field and the reflected field is determined by the distance between 
the dipole layer and the interface d and the radiation angle SI as 2/cidcosSl. Therefore the total radiation field at 
R(0, 0, R ) is the phase shifted sum of the E Rl and E R/2 fields in the EqGHand the EqCEQ We have 
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E x , total = -^—iki (n x cosn- u*sinfi - r 2p (u x cos 12 + sml2)e i2fcldcosa ) e l{klRo cosn -^) 
' iki (cosflL' xx fi x - smQL' zz iJ, z )e i( - klRoCOsn -^ 



a 2 Si 



■^ 7r ■]„ ( | _/fy_ e i2fcidcosf2 N j e i(kiRo cosfi-wt) 

v cos 51 s cos 12 / 



_ , Ztt ( 
E v , total = — — iki - 
a J £i V( 

a 1 e\ \ cos 12 y 

Ez ,total = ^ik, f_^ s inl2 + M2 ^+r 2s (^sinl2 + / , z ^)e^ dcos ^e^^ cosn -^ 
a^ei \ cos 12 cos 12 / 



27T - lkl (-sin 121^ + ) e ^«oco 8 o-. t ) ,,,, 



"'"i i » "-^xxf^x ' o 

a^ei \ cos 12 

Here the induced dipolc components /i x , [i y ,andii z are as defined in the Eq |12l and the L'^ are as defined in the 
Eqj5] Considering the fact that the thickness of the monolayer is generally much smaller than the magnitude of the 
optical wavelength, i.e. d <C A, the phase factor e ,2fcllicosn — > 1. Therefore, L' u becomes 

2nicosl2 t 2eifc 2z 

l^xx 



ri2 cos fli + n i cos Q t £2k\ z + £i&2z 
2nicosl2i 2k\ z 



«2z 

2ri2COsl2i 2e2fci 



m COS 12j + Tl2 cos 12* fci2 + &2 



712 cos 12^ + ni cos 12 t ^^lz+Elfe 



(36) 



These La values are clearly the results for the two phase model. This indicates explicitly that when dealing with the 
macroscopic electric field at the interface, the existence of the interfacial molecular monolayer has minimum influence 
on the macroscopic field. Therefore, the problem of the linear macroscopic dielectric constant for the interface layer 
is not a issue in this treatment. This shall be discussed in detail in the Section ITVl A. 

Now the radiation field components in the second harmonic frequency are 

j^2u) 27T 2u /„„„ o r 2w 2lu ■ r> r2w 2w\ tiki" Rn cos Q 2 iu —2ut) 

a z ei{2uj) 

t?2u> _ ^ 7r ,-i,2m f ^ t2u> 2u\ i(fef"fl cosQ 2M -2a;i) 

^ ~ a 2 ^)^ 1 \cos^ Lz ^y ) e 

El" = , 2 1 J kf (-BmSh u I% t % > + ^^L^n 2 A e <(^^»"..-^) (37) 
a z ei(2uj) \ cosS2 2 w ' / 



Here fc 2w = 2wy / ei(2a7)/ c, where c is the velocity of light in vacuum. Now with the second harmonic field components 
determined, the intensity of the SH field can be directly calculated by the magnitude of the time averaged Poynting 
vector as defined below 



I(u) = — v / iiM|£iM| (38) 

Then put the expressions in the Eq |3"T1 and the Eq[T2]mto the Eq[3Hl and after taking care of the proper projection 
coefficients of the fundamental and second harmonic electric fields, we have 

32tt 3 w 2 sec 2 12 2w , ,„ ., 2t2 , > /om 

1{2uj) = W^M IXe// :2w)l ( } (39) 

Xe// = pL(2w) • g(2w)] • x(2c^) : [L(w) • e»][L(u,) • e(u)} (40) 

Here is the intensity of the incident laser beam as defined in the Eql38l and 12 2w is outgoing angle of the 

second harmonic radiation from the surface normal. ei(w) and ei(2w) are the bulk dielectric constants of the upper 
bulk phase at the frequency to and 2u>, respectively. e(2u) and e(w) are the unit vectors of the electric field at 2uj 
and w, respectively. L(2oj) and L(w) are the tensorial Fresnel factors for 2u> and u>, respectively, as defined as in 
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the Eq[36] The scaler property Xe//(2w) is the effective macroscopic second-order susceptibility of the interface. 
The x(2o>) represents the area-averaged macroscopic second harmonic susceptibility tensors defined as Xijfc(2w) = 
lu(2u)ljj(u!)lkk(u)(3ijk/a* 2 = N s lu(2u)ljj((jj)lkk(u)([3i<jik<), with a as the lattice constant, N s as the surface number 
density, and ( ) represents ensemble average over different orientational distributions. 

The expressions in the Eq |37l and the Eq |39l are exactly the same as those in the previous treatment by Heinz 
and Shen.— ' 34 i 53 ' 54 i 55 This indicates that the treatment with the square point-dipole lattice model together with a 
molecular optics approach is at least phenomenologically correct. It also further validates the treatment in the 
original literatures. 

It can be shown accordingly that for the SFG from the interface, one has 

I(usf) = 8n u secsF | x eff (u S F)\ 2 I(vi)I(u) 2 ) (41) 

c j [ei(w S F)ei(wi)ei(w 2 )J 1/2 1 1 

Xeff = [M^sf) ■ e(uj S F)} ■ x(^sf) : [L(^i) ■ e(wi)][L(w 2 ) • e(w 2 )] (42) 
Now, the following general issues need to be discussed. 

a. Even though the expressions in the Eq l37l and the Eqj39]were obtained through the square induced dipole lattice 
model, using other type of plane lattice model gives the same asymptomatic expressions at the far field, except that 
the area-averaged macroscopic second harmonic susceptibility tensors take slightly different expressions. Certainly, 
the near field results have to be different for different lattice models, because generally no analytical asymptomatic 
results can be reached for the near field case. In other words, the molecular optics approach described here can be 
used to treat the near field cases, this is an clear advantage of the molecular optics approach. 

b. Here the area-averaged macroscopic second harmonic susceptibility tensors Xijkfad) an d the factor 1/a 2 in 
the Eq[37| needs to be discussed, a 2 is the average area per induced dipole in the lattice plane under the square 
induccd-dipole lattice model. It comes into the expression naturally because of the summation, or integration, over 
the whole lattice plane. Therefore, fii/a 2 is the area averaged induced dipole, and Xijk(2ui) is the area-averaged 
macroscopic second harmonic susceptibility tensors. In the previous treatments on second harmonic generational^ 
an infinitesimally thin polarization sheet layer was used, and the polarization for this interface layer was defined in the 
form of P(2uj) = x(2w)5(z) : E{uj)E(uj), or sometimes P(2u>) = x(2w) : E{oj)E(u)^^L In this way, it is a little 
confusing when one try to make the dimension analysis of the formulae when trying to apply the interface polarization 
term into the Maxwell equations. However, in this molecular optics treatment, there is no ambiguity involved, and 
the x(2w) tensor has to have the dimension of an area -averaged second order nonlinear polarizability. 

c. The ensemble average of the Xijfc(2w) tensors. In calculation of the summation and its asymptomatic results 
in the far field, the actual size of the lattice constant a can be treated arbitrarily as long as the lattice constant a 
is microscopic. Because the additive nature of the summation in the Eq |251 it is easy to show that even though the 
surface is not a square lattice for each individual dipole, one can always divide the surface layer into the square lattice 
with microscopic unit cells each containing a group of individual neighboring induced dipoles. Therefore, this approach 
actually allows making ensemble average with each unit cell. It even allow the unit cell to contain several layers of 
molecules close to the interface region. This further allows the inclusion of the contributions of the quadrupolar and 
the interfacial discontinuity terms into the effective unit induced dipole cells at the interface region i 4 i 26 ' 31 i 35 ' 97 i 98 Of 
course, each of these different treatments has to be involved with different microscopic description and symmetry 
analysis of the area-averaged effective microscopic polarizability tensor terms Xijk of the unit cell as defined. 

With above issues cleared, now we may say that the Eq[37] and the Eq(35] are general for the second harmonic 
generation from the whole interface region, with the possibility to make microscopic treatment or calculation with 
different microscopic models of the interface region. Similar treatment can be performed for the sum frequency 
generation process from the interface. Its detail is not discussed here. 

In the SHG, most of the time a rotationally isotropic molecular interface is studied. If we let a.i n be the polarization 
angle of the incident electric field, let 7 out be the polarization angle of the second harmonic electric field, Xef f,a in --y out 
be the corresponding effective macroscopic nonlinear susceptibility. Then the Xef f,a in ~^ out can be expressed into 
the linear combinations of the three independently measurable terms with the following polarization combinations: 
Si n — Pout ) 45°„ — s out , Pi n — Pout, with s as the polarization perpendicular to the incident plane, andp as the polarization 
in the incident plane. Let Cl u and f2 2w be the incoming and outgoing angles, respectively, one has 2 ^ 

Xe/i>„«-7o Ut = Xe//,45°- s sin7sin2a + x e //,s^pCOS7sin 2 Q! + Xe//,p-pCOS7Cos 2 Q! (43) 

Xeff,s- P = L zz (2cu)L yy (oj)L yy (uj)smn 
Xe//,45°-s = L yy (2uj)L zz (uj)Lyy(uj)sinn u jXyzy 
Xeff,p-p = ~2L xx (2uj)L zz (uj)L xx (uj) cosfl 2 Lj sin fl u cos fl^Xxzx 
+L zz (2u)L xx (uj)L xx (u) smn 2u j cos 2 Cl u Xzxx 

+L ZZ (2lo)L zz (uj)L zz (uj) sin 2w sin 2 ft^Xzzz (44) 
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Details of the quantitative experimental measurement and interpretation of the interface SHG, as well as the SFG- 
vibrational spectroscopy, of molecular interfaces with these formulations can be found in the recent literatures ! 21 ! 25 

IV. DISCUSSION 

Two of the central issues with general interests in the quantitative interpretation of the SHG and SFG-VS data 
are whether the molecular monolayer has a linear macroscopic dielectric constant and how the microscopic local 
field effects is evaluated. Actually previous researchers already provided most of the answers as reviewed in the 
introduction section above. Because the SHG and the SFG-VS have been proven as the sensitive probes for interfaces 
with the submonolayer cover age the treatment based on the more realistic discrete induced dipole lattice model 

is developed in this work. Based on the microscopic treatment of the SHG and SFG-VS from the molecular interface, 
these issues are to be discussed. 



A. The linear macroscopic dielectric constant of the molecular monolayer or submonolayer 

The concern of the linear macroscopic dielectric constant of the interface is from the beginning of the development 
of SHG and SFG-VS as the quantitative surface spectroscopic probes.— Here we quote the following sentences from 
the Chapter 'Surface Nonlinear Optics' in Shen's classical textbook) 2 ^ 

'In the usual simplified approach, the surface microscopic layer is assumed to have a characterized thickness with 
optical constants different from the bulk. Because of the linear transmission and reflection of the light at the boundary 
surface usually are dominated by the bulk properties, the surface layers have little effect on the linear wave propagation. 
As far as nonlinear optical effects are concerned, we can assume that the surface layer have the same linear refractive 
indices as the adjoining bulk media, but their nonlinear optical susceptibilities are different from the bulk. Unlike the 
linear case, the surface layer can strongly affect the nonlinear optical output in some cases.' 

The insight in these words is the basis for the later treatment by Shen to substitute the macroscopic dielectric 
constant of the interface layer in the early treatmen t 34 ' 54 with the microscopic local field factors. 41 ' 53 ' 56 This implies 
that the linear macroscopic dielectric constant of the interface layer does not need to be considered. Here the rigorous 
molecular optics treatment in the Section [H] and IIIII with the discrete induced dipole lattice model is providing a 
further support for this approach. 

Based on the following reasonings, the linear macroscopic dielectric constant of the monolayer or submonolayer is 
not an issue in the molecular optics approach described in the Section |TT] and IIIII explicitly. 

a. Firstly, we consider a few nonlinear induced dipolcs situated at the interface between the two isotropic bulk 
phases, for example, the vacuum/metal or other fluid/solid interfaces. These interface induced dipoles are situated in 
the vacuum phase and they does not even form an interface layer. Therefore, under the submonolayer condition up 
to a full monolayer condition, the treatment as presented in the Section [H] and IIIII with the discrete induced dipole 
lattice model is rigorous, and no linear macroscopic dielectric constant can be invoked for the molecular layer. When 
this is validated, then consider the case that there is a submonolayer up to a full monolayer on top of the first full 
monolayer at the air/metal interface. The treatment of this top layer must be the same, i.e., no need to invoke a linear 
macroscopic dielectric constant for this top layer, and only the microscopic local field effect, i.e. the dipole screening 
effect, need to be consideredi 56 ' 88 Since the interface layer in between any two isotropic bulk phases is usually one or 
few molecular monolayer thick, the discrete induced dipole lattice model must be a generally realistic treatment of the 
interface problem. In simple words, this implies that there is no such thing can be defined as the linear macroscopic 
dielectric constant for the interface layer. 

b. Secondly, let us consider the case for the submonolayer adsorption of organic molecules at the air/liquid, e.g. 
submonolayer of surfactant molecules at the air/ water interface. Since the air/water interface is known as sharp as 
less than 4A ) 99 i 100 i 1C11 one can surmise that the hydrophobic part is in the air phase and the hydrophilic part is in 
the water phase. In this case, the treatment in the Section [H] and IIIII with discrete induced dipole lattice model is 
generally a realistic treatment of the interface problem, especially for the SFG-VS when the tail and the head group of 
the intcrfacial surfactant molecule arc probed with different vibrational frequencies. The radiation from the molecular 
induced dipole in the phase 2 (£2) can also be treated similarly and readily as in the Section HT1 and Hill Nevertheless, 
in these cases, the continuous surface sheet layer model does not provide a realistic representation of the microscopic 
molecular picture. 

In the above cases, the linear macroscopic dielectric constant of the interface layer should not be a concern for the 
surface nonlinear optics whenever the surface contribution is dominant. Therefore, the macroscopic optical field at 
the interface can be well described with an explicit two phase model, as in the Eq |36l The major concern by Roy 57 
on the macroscopic linear anisotropy of the interface monolayer in the SHG treatment is not a real issue at least for 
the submonolayer up to the monolayer regime^ Nevertheless, according to the treatment presented in the Section 
HI1 and IIIII the local field effects of the induced dipoles at the interfaces need to be properly treated. This separation 
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of the macroscopic and the microscopic anisotropy of the interfacial monolayer is the ensurence of the simplicity and 
effectiveness in interpretation of the surface SHG and SFG-VS data. 

Of course, above discussion may not be fully suitable for the situation when the interface consists of a big number 
of layers of the anisotropic nonlinear induced dipoles. However, as long as such film is thin enough not to significantly 
alter the macroscopic reflection and transmission coefficients of light from the interface between the two isotropic 
phases, there is no reason that one should worry about the issue of the linear macroscopic dielectric constant of the 
film itself^&iSl 



B. The planewise dipole sum rule and the local field factor in the interface layer 

The treatment of the microscopic local field effect as expressed in the Eq[13]in the Section [ill is similar to the results 
as Ye and Shen's earlier workd The only difference is that now the general contribution from layers other than the 
interface monolayer is now included. To implement these formulae to correctly calculate or estimate the local field 
factors, the planewise dipole sum rule in the molecular crystal dielectric theory need to be invoked i 73 i 74 ' 75 i 76 ' 77 i 78 The 
implementation of the planewise dipole sum rule in studying the interface problem was discussed by Munn et al. in 
the 1990s i 49 i 79 ' 80 i 81 ' 82 However, its implication has not been picked up by others in the SHG and SFG-VS community 
so far. These issues needs to be addressed. 

Now let us look at the actual problems in the interpretation of the SHG measurement data. In their SHG studies 

of the self-assembled-monolayer (SAM) at the gold substrate, Eisert et al. carefully compared the results with the 

macroscopic three layer model and the two layer model, as well as the local field corrections. They concluded that 

using the two-layer model without local-field correction gave most satisfactory agreement of the molecular orientation 

with the results from the NEXAFS and IR spectroscopy measurements^ According to Eisert et al., the parameters 

used to calculate the local field factors in the Eq[T3] and the EqQ3] of the pNA [O2N-C6H4-] monolayer in the SHG 

3 3 

experiment are: a = 5 to 5.5A, a aa = abb = 16. 4A , a cc = 27.3A , and 8 ~ 52° from NEXAFS measurement. Then, 

□ 3 o 3 

one has a xx = a yy = 19. 8A and a zz — 20. 5A , then l xx — l yy — 3.51 to 2.16 and l zz = 0.403 to 0.473. Therefore, 
the effective microscopic dielectric constant defined as e' — n' 2 — l yy /l zz becomes 8.71 to 4.57, or n' = 2.95 to 2.14. 
Here the definition e' — n' 2 — l xx /l zz follows the definition by Wei et al.. 41 These are unreasonably large values and 
the calculation by Eisert et al. concluded that in order to get the ~ 52° from the SHG data, n' has to be close to 
unity, i.e. neglecting the microscopic local field effect. 

Such conclusion by Eisert et al. is in direct disagreement with the formulation provided by Shen et al^ How- 
ever, the formulation with the microscopic local field factors by Shen et al. is validated by the treatment in this 
work. The answer to Eisert et aZ.'s difficulties in calculating the local field factors lies in the planewise dipole sum 
rule ) 73 ' 74 i 75 ' 76 i 77 ' 78 which Munn et al. used to tackle the problem of overestimation of the local field factors by putting 
the linear polarizability of the whole molecule into the Eq[T3 l 49 ' 79 i 80 ' 81 i 82 

The planewise dipole sum rule was originally developed in explaining optical and dielectric phenomena in molecular 
crystals, such as the long range coupling of the exciton in the molecular crystal, in the early 1970si 73 i 74 ' 75 i 76 ' 77 ' 78 It 
has been known that the plane sum in the summation of the dipolc-dipole interaction as in the Eq(6j which leads 
to expressions of the local field factors in the EqfTSI falls off exponentially as the perpendicular distance from the 
plane of the origin increases^ 7 ' 7 Philpott et al. showed that even for the strong dipoles in the molecular crystals, the 
contribution of the immediate neighboring layer is only less than 1% of the contribution of the plane of the origin, 
and the total sum of all the rest of the layers up to infinity distance other than the plan of origin is generally less 
than 10% 

Based on this planewise dipole sum rule, when the lateral distance between the molecules is smaller than the 
length of the molecule, there is no reason to view the molecule as a whole when calculating the local field factors 
using the Eq |13l Munn et al. thus proposed a bead model to divide the long chain into a chain of sphere beads in 
calculating the local field factor for each submolecular segment layer for the long chain molecules in the Langmuir- 
Blodgett fihnj 49 i 79 i 80 i 81 i 82 This bead model is also unknowingly supported by the calculations by Ye and Shen on the 
contribution of the local field factors from the image induced dipoles at a vacuum/metal interfaced They conclude 
that for Emetai = 10, i.e. {e — l)/(e + 1) ~ 1, when d/a > 1/2, the image induced dipole contribution is negligible. 
Here d/a > 1/2 is just the sphere bead condition in Munn's treatment if the interfacial induced dipole and its image 
induced dipole are considered as a whole unit. 

In our treatment, the planewise dipole sum rule is explicitly described with the Eq llOl Evaluation of the £q term 
indicates that in most of the cases, even the immediate neighboring term, i.e. / = ±1, is often negligible when d is 
larger than 3A. This indicates that when dealing with the local field factor calculations with the EqQ21 the depth of 
the segmentation of the molecular layer is generally not more than 3A. However, this number may not be correct for 
the big chromophores with conjugate structures, which may require quantum mechanical treatment or more complex 
models rather than the simple classical dipole model. 

With the planewise dipole sum rule in mind, now one can understand why in the Eisert et al. case^i that the local 
field factors are over estimated than the actual value. Because the SAM monolayer is closely packed, if one use the 
an value of the whole pNA group, the la factors should be more deviated from the unity value. If one use the sphere 
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bead model, the la is going to be much more close to unity. Actually, if closely inspect the formulae in the EqfT3l 
one can expect that when | aa^o/2a 3 \> 1, the l xx and l yy can even become a large negative value. If the sphere 
bead model is not used, a closely packed long chain molecular monolayer would easily have unreasonable local field 
factors according to the Eq |13l without considering the planewise sum rule. Therefore, the planewise dipole sum rule 
is very important and has to be implemented for the evaluation of the microscopic local files factors in the molecular 
monolayer. 

One direct conclusion from the segmentation of the interfacial molecules for the implementation of this planewise 
dipole sum rule is that because the different groups in the same molecule has different linear polarizabilities, and 
because the same molecular group may have different polarizabilities at different frequencies, their local field factors 
can be significantly different when these groups occupies different segmented layers when the monolayer is closely 
packed. This immediately poses questions to the practice of using the same and simple value for the interface local 
field factors in the SHG and SFG-VS studies^^^ One of the example is for the SHG and SFG-VS studies of the 5CT 
[CH3-(Cti2)4-(C6H4)3-CN] monolayer at the air/water interfaced In the interpretation of the SHG measurement of 
the -(CeELi^-CN chromophore and the SFG-VS measurement of the -CH3 and -CN groups of the 5CT monolayer, the 
same local field factor value, i.e. n' = 1.18, was used for all three molecular chromophore or groups. Close inspection 
shows that the orientational angle thus obtained is actually inconsistent with other studies. 102 ' 103 ' 105 ' 106 Therefore, 
whether the same n' = 1.18 value is suitable for the three groups needs to be re-examined. In a recent study in 
our research group, using a self-consistent approach we found that the n'(800nm) ~ 1.5 and n'(400nm) ~ 2.5 for 
the -(CgEU^-CN chromophore from the SHG data, significantly different from that of the -CH 3 or -CN group, and 
consistent with the large polarizability of the -(CeH4)3-CN chromophore i 104 ' 105 The orientational angle thus obtained 
for the -(C6H4) 3 -CN chromophore in the closely packed Langmuir monolayer at the air/water interface is about 20°, 
consistent with the orientational angle obtained in other studies i 102 ' 103 

Now, if the e' or n' values should be different, why in the past SFG-VS studies it was quite successful to use the 
simple values such as n' = 1.18, or the average value of the optical constant of the two neighboring bulk phases for 
molecular groups and chromophores at the common dielectric interfaces ? 21 ' 25 ' 53 ' 58 ' 68 ' 106 Here we can show that these 
values are actually quite reasonable approximation of the common simple molecular groups such as -CH3, -C=0, etc. 
at the interface. 

The experimental linear polarizability values of many common molecular group are as listed in the Table [J 93 ' 94 ' 95 

3 

The values for simple organic groups are typically in the range of 1.5 to 3.0 A . Considering that the segmentation 

depth is expected to be not more than 3A for simple molecules with no big chromophores with conjugated structure, 

° 3 

as a realistic approximation, the typical a value used in the calculation is set in between 2 to 4 A . Considering 

o 2 

the fact that at the molecular interface a typical molecular group per area is about 20 A (square lattice constant 

a ~ 4.5A) to 30 A 2 (o ~ 5.5A), simulations in Table HT1 shown that the typical e' value for these groups are 1.17-1.56, 
i.e. n'=1.08-1.32, depending on the lattice constant and the value of the polarizability. 

In above calculation we assumed that the monolayer is closely packed. If the interfacial monolayer is sparse, the 
local field effect shall be small. For the situation of submonolayer adsorption at the liquid-liquid interface, then the 
dipole interactions of the solvent molecules in the interface plane need to be considered. 

These values are generally in between the vacuum and the substrate dielectric constant, as used in the literature. 
They are also close to the estimation made by Zhuang et al. using a modified Lorentz model of the interface;^ 3 - and 
the simple average of the optical constant of the two adjacent isotropic bulk phases! 52 ' 58 ' 68 However, even though 
these values worked almost fine, the models and the reasonings behind it prevented the treatment of the molecular 
details of the molecular interface. With the development of the laser techniques and the measurement methodology 
in the surface SHG and SFG-VS ) 21 ' 25 ' 39 ' 40 now such details can be interrogated from the careful experiments. This 
shall provide new opportunities for the interface studies. 



TABLE I: Experimental polarizabilities values for common molecular groups from the literatures.— 



group 


a(X") 


group 


a (A a ) 


-COOH 


2.86 


-SH 


3.47 


-CH 3 


2.24 


-CN 


2.16 


-CH 2 - 


1.84 


-NH 2 


1.76 


-C=0 


1.82 


H 2 


1.45 



Of course the above simple estimation does not include the bigger chromophore, such as phenyl, and biphenyl etc. 
Their polarizability is generally much larger, and the local field effects can be much stronger when they are closely 
aligned at the interface. According to the Eq |131 it is possible to have the l xx = l yv to be smaller than unity or even 
go to negative values, and l zz to be larger than unity, when the molecules with relatively larger linear polarizability 
are closely packed even far from resonance. These are the unique phenomena for the two dimensional anisotropic 
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TABLE II: Simulation of local field factors e and n with the a value between 2 and 4 A 3 . The case for a = 1.5A 3 is the 
simulation results of the interfacial water molecules. 

° 

ct(A ) a(A) l xx = l yy l zz £_ n' 

2.0 4.5 1.11 0.84 1.33 1.15 

2.0 5.5 1.06 0.90 1.17 1.08 

4.0 4.5 1.25 0.72 1.74 1.32 

4.0 5.5 1.12 0.82 1.37 1.17 

1.5 3.5 1.19 0.76 1.56 1.25 

1.5 4.0 1.12 0.83 1.36 1.16 



local field factors. According to the Lorentz relation, this is not possible for the optical bulk materials with normal 
dispersion. Such phenomena are certainly yet to be explored experimentally. 

The local field factors in the molecular layer is dependent on the orientation and also on how the segmentation is 
done, and one generally does not have the detailed knowledge of these information beforehand. Therefore, different 
approximation approaches were employed in the past SHG and SFG-VS studies. Similar to what Munn et al. did 
previously , 49 ' 79 i 80 ' 81 i 82 we invoked the planewise dipole sum rule here to show that in calculating the local filed factors 
using the microscopic point dipole model using formulae in the Eqll3[ extreme care must be taken in order not 
to overestimate the local field effects. However, this does not imply that the sphere bead model always generates 
more reasonable results for the local field factors. We want to point out that when the square lattice constant is 
close to or even smaller than the size of the molecule under studying, the simple point dipole model is expect to 
breakdown, and the planewise dipole sum rule can be used as the remedy to some extent. One expects that there is 
no simple rule to segment the chain molecules in calculating the local field factors. Down to the detailed molecular 
level, quantum mechanic treatment of the electron density and their polarizability should be employed to make more 
accurate description of the molecules at the interface. Nevertheless, the above estimation ca give the reasonable upper 
and lower bound of the local field factors. 



C. SHG or SFG-VS experiments to test the microscopic model 



Here we discuss how SHG and SFG-VS experiment can be used to test the microscopic model. 

The key point raised in this study is that the molecular optics treatment of the interface nonlinear radiation can 
provide a detailed microscopic description of the nonlinear optical processes and the molecular behaviors at the 
interface. The key to test the microscopic model is to find way to determine the microscopic local field factors or to 
quantitatively determine their effects in the SHG or SFG-VS experiment, and to compare these experimental effects 
with the molecular details in the molecular monolayer. 

There were few attempts trying to experimentally determine the interfacial effective optical constants with SH C 58 i 68 
and SFG-VS. 67 All of them concluded for a value close to the simple arithmetic average of the optical constant of the 
two adjacent bulk phases. These works were not based on the microscopic model and they can not be used to test 
the microscopic model. Because the microscopic model predicts that the local field factors, which is connected to the 
effective optical constant of the interface layer through the relationship e' = n' 2 = lyy/hz i 41 ^ 3 must depend on the 
detailed structure of the molecular monolayer and may have values other than the simple arithmetic average or the 
prediction with the modified macroscopic Lorentz model of the interfaced 

The first test is to determine the local field factor l yy /l zz for the closely packed Langmuir monolayers at the 

° 3 o 2 

air/water interface with large chromophores. According to the EqfTBl, when an ~ 15A and a 2 ~ AO A , one has 
that e' = lyy/lzz ~ 2.10, thus n' ~ 1.45. These kind of conditions can be easily satisfied with chromophore with 
two or three phenyl groups, such as the alkyl cyano-biphenyl (nCB) or the alkyl cyano-triphenyl (nCT) Langmuir 
films. It has been shown that in the closely packed 8CB Langmuir monolayer at the air/water interface, since the 
molecular surface density is always known, a self-consistent analysis of the polarization dependent SHG data at 
different surface density can give the molecular orientation angle and the e' value at different surface densities ! 37 ' 105 
Such self-consistent analysis was suggested by Munn et al. before, because they realized that the determination of the 
molecular orientational angle and the determination of the interface local field factor are interconnected^ Such self- 
consistent analysis can only be carried out for the Langmuir or Langmuir-Blodgett monolayers because their surface 
density are known and can be varied controllably. It turned out in our SHG studies that for the 8CB monolayer, the e' 
changes from 1.7 to 2.4 when the surface density changes from 51 A to 39 A per moleculei 104 i 105 These results agree 
very well with the EqfLU and the orientational angle thus obtained agree well with the rod model of the chromophore. 
This is a very convincing example to show that the local field factors of the molecular monolayer do not have simple 
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values close to the arithmetic average of the two adjacent bulk phases. The detail of this study is to be published 
elsewhere, 10 V 05 

Another test can be designed with the SFG-VS measurement of the small rigid linear molecules, such as the 
acetonitrile molecule, at the air/liquid interface. For example, SFG-VS can selectively measure polarization dependent 
vibrational spectra of both the —CH3 and the —CN groups of the interface acetonitrile molecule* 1 - 7 Since the molecule 
is linear, the orientation angles, which are dependent on the values of s'_ CN and s'_ CH3 respectively, determined from 
the —CH3 and the —CN SFG-VS data should give the same value. Thus the ratio e'^cn I e '-CH 3 can ^ e determined 
experimentally. Using the polarization null angle method in SFG-VS ) 21 i 108 i 109 i 110 the ratio of the local field factor for 
the —CH3 and the —CN groups, i.e. s'_ CN j e'_ CH ^ can be obtained quite accurately. This value is indeed different 
from unity and it is beyond the experimental error bar. This observation can not be explained by the previous models, 
but can be well explained and quantitatively analyzed with the microscopic model in this work. The detail of the 
experimental study and analysis is going to be reported elsewhere* 1 ^ 1 - 

More SHG or SFG-VS experiments can be designed accordingly to test the microscopic model. Detail analysis 
of SHG and SFG-VS data is now possible with the recent development of quantitative polarization measurement in 
SHG and SFG-VS i 21 i 25 ' 112 ' 113 These studies not only can provide test for the microscopic model, but also can provide 
detailed information of the molecular interface according to the microscopic model. 



V. CONCLUSION 



Because the SHG and the SFG-VS have been proven as the sensitive probes for interfaces with the submonolayer 
coverage^ 4 ^^ the treatment based on the more realistic discrete induced dipole lattice model is developed in this 
report. Recent development in quantitative polarization and symmetry analysis in SHG and SFG-VS have provided 
new opportunities in the interface studies, 21 ! 25 These development also provided more accurate data to test the SHG 
and SFG-VS theories at the detailed molecular level -S^^i 5 ^ 8 . 

In this report we have taken the molecular optics theory approach and treated the molecular interface with more 
realistic discrete induced dipole lattice model. Based on this model, the following results are derived: a. the detailed 
expressions of the local field factors in the interface molecular layer; b. the detailed expressions for the far field radiation 
of the SHG as well as the SFG-VS process from the interface induced dipoles. It turned out that the asymptotic results 
for the far field radiation is in the same form as the results derived from the famous infinitesimally thin polarization 
sheet layer by Heinz and Shen using the Maxwell equation with the boundary conditions* 4 ^ as well as their later 
modifications with the consideration of the microscopic local field factors using the classical dipole mode h 41 ' 53 i 56 
This not only validates the effectiveness of the original formulation from the works by Shen and colleagues, but also 
validates the success of the molecular optics theory approach developed here. 

According to Born and Wolf, the molecular optics theory can directly connect the macroscopic optical phenomena 
to the molecular properties, and can provide deeper physical insight into electromagnetic interaction problems than 
does the rather formal approach based on Maxwell's phenomenological equations, 83 ! 84 Based on the microscopic and 
discrete induced dipole lattice model, the problem of the macroscopic dielectric constants of the interfacial molecular 
monolayer is discussed and clarified. We explicitly demonstrated that the macroscopic dielectric constant of the 
molecular submonolayer and monolayer does not need to be invoked in the microscopic model of the surface nonlinear 
optics. Based on the planewise dipole sum rule and the previous work by Munn et ai, the issue on how the microscopic 
local field factors are evaluated was discussed. According to these results, the effectiveness and limit of the simple 
models on the effective dielectric constant of the interfacial molecular monolayer arc discussed and evaluated. However, 
for simple and small molecular groups, the simple models does provide good approximation of the microscopic local 
field effects of the interfacial molecular monolayer. SHG and SFG-VS experimental tests of the microscopic discrete 
induced dipole lattice model are also discussed. With these developments, many previous SHG and SFG-VS studies 
are better understood and evaluated. Moreover, these developments can provide a full microscopic description of the 
nonlinear radiation from the molecular interface. Detailed molecular information can be obtained from the model 
developed in this work, together with better and more accurate polarization measurement data in SHG and SFG-VS. 

The microscopic molecular optics approach can also be applied to the treatment of the ellipsometry response 
from the molecular monolayer interface. When ellipsometry is to be applied to study such interfaces, a microscopic 
molecular theory for the ellipsometry is certainly needed. So far, the treatment on the ellipsometry and the treatment 
of the SHG or SFG-VS are not fully consistent in description of the anisotropy in the molecular layer. We shall discuss 
these issues elsewhere. 

In conclusion, aside from the macroscopic theory of the surface SHG and SFG-VS, an effective microscopic molecular 
optics theory is developed in this work. Such development may shed new light on both the linear and the nonlinear 
optics interface studies. 
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